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SUMMARY 


This report describes a nonintrusive nuclear technique for monitoring fluid 
contents in sealed vessels, regardless of the fluid distribution inside the 
vessels. The technique is applicable to all-g environments. It is based on the 
differences in Cesium-137 gamma ray attenuation coefficients in air and the test 
liquids . 


INTRODUCTION 

We recently developed a nuclear technique for monitoring water levels 
in pressurized vessels. It was based on the differences in attenuation coefficients 
of water and air for Cs 137 (662 Kev) gamma rays. It is now proposed that the 
nuclear monitor can be adapted for measuring the total fluid content in closed 
vessels regardless of the fluid location inside the vessels. Such a monitor would 
be useful for 0-g environment or other situations where effective accelerations 
could be different from that on earth. However, such a monitor would require 
specialized source and detector geometries. 


In this paper, we describe a system for monitoring the arbitrarily distributed 
liquid content of a vessel and theoretically demonstrate its effectiveness using 
a computer simulation. 


CONFIGURATIONAL DETAILS 

The principle of operation of the standard nuclear monitor has been described 
in detail in reference 1. Briefly, it depends on the fact that 


I (Air) » I (Liquid) 


where I (Air) = I e X ° + Xa ^ 

t o 


and 1^ (Liquid) = I e 


- ( 2 U C * c + V*> 


where I = Incident Gamma Ray Intensity 

= Linear Attenuation Coefficient for the Container Material 

M = Linear Attenuation Coefficient for Air 
a 

= Linear Attenuation Coefficient for the Test Liquid 

x = Container Wall Thickness 
c 


( 1 ) 

( 2 ) 







x q = Gamma Ray Pathlength in Air 
= Gamma Ray Pathlength in Liquid 


The system in reference 1 utilizes a collimated source of gamma rays measured by a 
detector having wide acceptance angle. This approach is suitable for the more 
conventional applications where the fluid seeks a stabilized position in the vessel. 


This paper considers a redesigned system such that the monitor would work 
regardless of the fluid's location/distribution inside the container. Several 

source and detector configurations were considered ( '. Some of them are listed 
below: 


1. Uniformly distributed, highly collimated point sources and detectors on 
opposing faces of the container. 

, . , 2 * U nif °rm linear source on one face with point detectors uniformly distri- 

buted on the opposing face of the container. 

3. Non-uniform linear source on one face with a uniform planar detector on 
the opposite face. 

4 * Uniform linear source, with various degrees of collimation.on one face 
with a uniform planar detector on the opposite face of the container. 

5. Uniform planar source on one face with a uniform planar detector on the 
opposite face of the container. 

. Aft ® r trying various source/detector arrangements, two major modifications 
to the scheme in reference 1 were made: 

1. The well-collimated point radiation source was replaced by a planar 
source uniformly deposited on one side of the container. (It was observed that 
the source could just as well be in the form of a narrow linear strip passing 
vertically through the center of one side of the container, if the fluid were to 
be distributed randomly inside the container. However, if the fluid were pressed 
against one of the walls or confined to corner regions due to centrifugal forces, 
l would be preferrable to have a planar source covering the entire side.) 

. , . 2 * , The s f cond modification required replacing the finite-size cylindrical 

detector by a planar detector covering the entire opposite wall of the container. 

Hpi-p-m Fl ® ure 1 ®^ows a schematic diagram of the proposed concept. The computational 
details for testing the system are described in the following section. 


COMPUTATIONAL PROGRAM 

The computer program CRANK is written in FORTRAN IV language for the 
Control Data CYBER 170 series digita l computer system with Network Operating 

(*) ' “ ‘ 

For the sake of specificity, we assume the container to be a SL x 2d x 2d 
rectangular box. 
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System (NOS) 1.4. The program requires approximately 52,000 octal locations of 
core storage . A typical case in which 40,000 trajectories are computed requires 
approximately 470 CPU seconds on CYBER 175. 

The program models the stochastic process in which a gamma ray traverses a 
medium subject to multiple Compton scattering(2) . The technique used for tracking 
the progress of each gamma ray through the medium has been reported in the litera- 
ture(3,4) and described in detail in reference 1* The medium modeled by the 
. program, is representative of a closed container holding variable amounts of a fluid 
located randomly within the container. 

The interior of the container is illustrated in figure 2. The region,- of 
exterior dimensions l x 2d x 2d, is subdivided into smaller volumes such that 
the cross section in any plane parallel to the x-z plane is an n x n uniform 
grid. The parameters S,, d, and n are defined by the user as program input. Within 
this • grid , the program assigns to, and identifies each cell with, an integer index 
between 1 and the total number of cells. If a cell is designated to contain fluid, 
the electron density (and thus the scattering properties) for that cell is assigned 
a value corresponding to the fluid. All cells are initially assigned an electron 
density corresponding to air, 3.65 x 1020 electrons/cm 3 . 

As currently configured, the program is subject to several limitations, each 
of which can be changed with minor modifications to the code. The shape of the 
container is assumed to be a parallelopiped with a square cross section in the 
x-y plane. The initial position of each gamma ray is assumed to be uniformly 
distributed in the x-y plane and to be contained within the container interior 
(i.e. | x | '< d, | y | « d) . Each gamma ray crossing the z = i plane within the 
container interior is assumed to be counted. 

Program execution is performed in sets where the number of sets is provided 
by the user. A single test contains repetitions for several fluid distributions 
where the number of distributions is a function of the grid selected by the user. 

For a 10 x 10 grid, 11 distributions are considered corresponding to fluid content 
of 0, 10, 20, ..., 90, 100/. For a 5 x 5 grid, 6 distributions are considered 
corresponding to fluid content of 0, 20, 40, 60, 80, 100%. For each distribution 
the number of gamma ray paths to be traced is specified by the user. 

-For each gamma ray path to be traced the program first assigns an initial 
gamma ray position and direction. The initial Z coordinate of the particle (ZI) is 
set to zero. The X and Y coordinates (XI, YI) are assigned uniform random values 
between -d and d. The initial direction of the particle is defined by the standard 
spherical coordinates 0 and <f) where the initial value of 0 is assigned a 
uniform random value between -0 L and 0 L and where 0 l is a limiting angle pro- 
vided by the user' ). The initial value of <f> is assigned a value randomly selected 
between 0 and 2ir. Following the technique of reference 1, the trajectory inter- 
section with each sampling plane is determined and the coordinates of this inter- 
section are found to either be contained in a cell containing fluid or an empty cell. 
The scattering at this point is thus related to the cell contents. If the Z coordi- 
nate of the particle becomes negative or if the the magnitude of the x or y coordinates 


W) 


e L could of course be as high as ir/2. However, calculations indicated that no 
significant advantage was gained for 0 > 30°. 
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exceeds d, the particle is dropped from further consideration. Once the Z 
coordinate of the particle exceeds 1 , the particle is counted. This process is 
repeated for each trajectory in each configuration in each test. 


Program Input 

All program input is accomplished using FORTRAN list directed reads. The 
data may appear anywhere in the field, and when more than one item is specified, 
the data items are separated by commas. The first items in a data set are the 
initial gamma ray energy and the cutoff energy in keV. If, as a result of scatter- 
ing, the gamma ray energy drops below the cutoff value, it is not counted by the 
detector. The next data items are the dimension parameters of the container 
interior, S, and d, in cm. The next data items are the number of sets, the number 
of grid lines (5 or 10) , and the number of trajectories per distribution. These 
data are followed by the electron density of the fluid in electron/cm^, then the 
distance between sampling planes in cm, and finally the limiting value of 9 for 
particle initialization, in degrees. If additional cases are desired in a single 
program execution, the entire data set must be repeated for each case. Sample 
data for a single execution to analyze two cases ( / are illustrated below: 

662., 200. 

30.48, 10.16 
10, 5, 400 
3 . 346E+23 
1.016 
30. 


662. , 200. 
30.48, 10.16 
10, 10, 400 
3.346E+23 
1.016 
30. 


In addition to the data described above, the program also requires a table 
containing probability versus scattering angle (0), as described in reference 1. 
These data should exist as an unformatted file on unit TAPE2. 


Program Output 

The program generates formatted output on two files. File TAPE6 contains 
a record of the program input followed by a graphic display of the fluid distri- 
bution and the transmission probability for each distribution. An example of this 


(*) 

For illustrative purposes, the fluid is assumed to be water. The container 
dimensions are taken to be 12" x 8" x 8" and the grids considered are 
(5 x 5) and (10 x 10). 


4 



output is contained in Table I,, for a 10 x 10 grid. Each distribution from no 
fluid to 100^ fluid is illustrated for each set, with the transmission for each 
distribution printed below the corresponding display. In the display, each "X" 
represents a cell containing fluid with embedded blanks representing cells contain- 
mg air. Following this display for each set, a summary over all sets is output 
\ s " icated in Table IT. This summary contains the average transmission and its 
{Standard deviation for each distribution. File TAPE1 provides a brief summary for 
each distribution including the number of counts, the number not counted, and the 
num er of particles dropped for x, y, or z coordinates outside the region or for 
energy below the cutoff value of 200 keV. 

In addition to the formatted output, the program generates graphic output 
presenting the transmission versus fluid content. An example of this output is 
included as figure 3. The solid line in this figure shows the average transmission 
over all sets and the dots show the results for each set. 


Program Listing 

A listing of the program and required subroutines follows. For completeness 
the program for generating gamma ray scattering probability vs. 0 is also included, 
reas of code reflecting limitations of the current program are highlighted by 

ConmiPnf-Q o cd j 


Pc npp A A* Tr>i$T?C r>|jTPl*T , T*'P t ! T, TA PFf "OUTPUT, TA P E I * 3 5«IKPU T * 

1 TiPr?) 

rrKMnK! PI, PC, At PMA 
FVTFPKAL Ft'*! 

PI « AC0S (-1 . ) 
p 0 * 2.P1PP-13 
At pma = 

TSCAT = GLFC1«(0.,pt*FUK’) 

W p TTP(f,1 ) 

1 FHPf'AT nMl,10Y,5HC7FRn>10y*5HALPHA» 

1 ?v,?7HT n T SCAT CROSS SFCTJ0N CPT) ///) 

ITF ( f , ? ) PC, ALPHA, TSCAT 
? FpP^ATnw0,lPF15.4,lPF??.4,lPF?9.A) 

NI.'M * ] PC, 

no ion i=i,ku m 

TH c TA = PI*FL HAT ( 1-1 )/l PC, 

THFTAP » THFTA + PI/3PC. 

PSCAT - GLFG15(TMFTA>TMFTAP,FUN) 

VP I TF ( ? 3 T p f t a , PSCAT 
VPITF(A,a) THFTA, PSCAT 

3 P 0 R H A T C?F1f,P) 

ICO CCMTIH’F 
CnMTTki'F 
STOP 

FA'O 
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no r> r> o r> 


^PROGRAM CRANK ( OUTPUT# IN PUT # TAPE 5 ■INPUT# TAPE6»0UTPUT# TAP E2# 

COMMON PI#RO# ALPHA# EZER0#EFIN4L#CLEN,C0I AM,DI ST,THETA L, XNZW 
DIMENSION TTT ( 100 ) L ' ^ 

DIMENSION XX(180)#YY(180) .i 

DIMENSION TRANS (100# 10) l 

k C0MM0N/RSLT/AVG(11)#STD(11) E 

f DIMENSION 16(100) 

1 F0RMAT(1H1#3X»*EZER0#EFINAL#CLEN#CRAD»NS5T#NGRD#NGAM# DIST#THETAL* ' 
1 *#XNZW*/ME16.8//3I10//3E16.8) M L 

2 i F ^ M lIIJ H1 nl X /*^!! ,f1ARV * ///lX,6X, * PERCENT WATER* »3X,*TRANSMISS ION*# 

1 oXj*bTD# DEV#*//) 

3 F0RMAT<5X#3F15.4) 

CALL PSEUDO 

PI • ACOSC-l.) 

GET SCATTERING DISTRIBUTION 

DO 10 1*1# 180 
READ (2) XX(I)#YY(I) 

10 CONTINUE 
SUM • 0. 

GET CUMULATIVE DISTRIBUTION 

DO 20 I*l#180 
SUM • SUM ♦ YY ( I ) 

YY(I) « SUM 

XX(I) • XX (I > ♦ PI/360. 

20 CONTINUE 
C 

C NORMALIZE IT 

C 

DO 30 1*1# 1B0 
’ VY(I) . YY( I ) /SUM 
30 CONTINUE 

RO « 2.818E-13 
C 

r READ INITIAL AND CUTOFF ENERGIES 

C 

40 READ (5#* ) FZERO# EFINAL 
IP(EOF( 5I.NE.0) GO TO 999 
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c 

c 

c 

c 

c 


* 


r 

* c 
c 
c 
c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 


c 

c 

c 

c 

c 

c 


READ LENGTH AND ‘RADIUS* OF CONTAINER 
FOR A RECTANGULAR CONTAINER THE 
•RADIUS* IS HALF THE WIOTH 

READ ( 5# * ) CLEN* CRAD 
CDIAM ■ 2.*CRAD 

READ THE NUMBER OF SETS# THE NUMBER 
OF GRIDS* AND THE NUMBER OF GAMMAS 
PER SET 

A 

READ ( 5# * ) NSET # NGRD# NGAM 

READ THE ELECTRON DENSITY OF 
THE FLUID 

READ<5*+) XNZW 

RFAD THE DISTANCE BETWEEN 
SAMPLING PLANES 

READ { 5# * ) DIST 

READ THE COLLIMATION ANGLE 

R FAD ( 5# * ) THETAL 

WRITE (6*1) EZERO#EFINAL*CL EN#CRAD#NSET*NGRD*NGAM#OIST*THETAL 
1 * XNZ W 

THFTAL - THETAL*PI/180. 

ALPHA - EZERO/511. 

CALL CALPLT(1.5#1.5#-3) 

NGSO • NGR D*NGRD 
NGRDP - NGRD ♦ 1 
NGRDM ■ NGRD - 1 

THE 100 LOOP IS OVER SFTS 

DO 100 K«1#NSET 

ZERO THE GRID 
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DO 50 I • 1* NGS 0 
50 IGC I > « 0 
C 

C THE 80 LOOP IS OVER FLUID CONCENTRATIONS 

C 

DO 80 I»1#NGRDM 
C 

C RANDOMLY FILL PART OF THE CONTAINER 

C 

DO 70 J-1#NGRD 
60 R « RANF(T) 

NP » NGSQ+R + 0.01 
IF (NR.LE.O) NR • 1 
IF(NR.GT.NGSO) NR - NGSO 
TFCIG(NR).EQ.l) GO TO 60 
IG(NR) « 1 
70 CONTINUE 
C 

C PROCESS THIS CONCENTRATION 

C 

CALL CQNFIG(IG# I# TTT# TR ANS#NSET# NGA M# NGRD# XX# YY ) 

80 CONTINUE 

00 go I»1#NGSQ 
IG ( I ) • 1 
00 CONTINUE 
C 

C PROCESS THE FILLED CONTAINER 

C 

CALL CONFIG dG# NGRD#TTT#TRANS#NSET# NGAM»NGRD» XX# YY ) 
100 CONTINUE 
r 

C DISPLAY THE RESULTS 

r 

CALL 01 SP ( TTT # TRANS* NS ET#NGRD) 

C 

C PRINT THE SUMMARY TABLE 

C 

WPITE(6#2) 

DO 110 I»1#NGRDP 

PER ■ (100./FL0AT(NGRD))+FL0AT(I— 1) 

WR I TE (6# 3 ) PER# AVG( I ) # S TDC I ) 

110 CONTINUE 
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GO TO 40 
909 CONTINUE 

CALL C ALPLT (0*> 0»> 999) 

STOP 

END 



c 

c 

c 

c 


c 

c 

c 

c 


c 

c 

c 

r 

c 

c 

c 


SUBROUTINE CONF IG( 16, IND#TTT# TRANS, NSET# NGAM,NGRD# XX# YY I 

CONFIG IS ENTERED ONCE FOR EACH 
FLUID CONCENTRATION FDR EACH SET 


COMMON PI 
DIMENSION 
DIMENSION 
DIMENSION 
DIMENSION 
DIMENSION 


#RO# ALPHA# EZERO#E FINAL# CLE 
XX( 180)#YY (180) 

TTT(IOO) 

NZERO(10)#NONE(10),TT(10) 

IG(NGRD#NGRD)#KG(100#10) 

LINE(IOO) 


N#CDIAM,DIST#THETAL#XNZW 


DIMENSION IZERO(10)#TRANS(100#10) 
EXTERNAL FUN 
DATA ISET/O/ 


DATA IZER0/10nH / 

1 FORMATdX# UdOAIflX)) 

2 F0RMAT(1H0#11(10A1#1X)) 

3 FORMAT! 1H0# 3X#*I#NHIT#NMISS#NNX# NNY#NNZ#NNE*// 
1 7110//) 

4 F0RMAT<1H0#F9.6#2X#10<F9.6#2X)) 

CP AD • CDIAM/2. 


STORE CONCENTRATION DATA FOR 
LATER DISPLAY 


NGRDP - NGRD ♦ 1 

NGRDM » NGRD - 1 

NFS 0 « NGRD*NGRD 

DO 10 1*1# NGRD 

JND ■ NGROMIND - 1) + I 

DO 10 J«1#NGRD 

KG(JND»J) - I G ( I# J ) 

10 CONTINUE 


RETURN IF NOT FULL 

IF ( IND.NE.NGRD ) GO TO IPO 

LAST CONCENTRATION (FULL)# 
DISPLAY AND PROCESS 

I SFT - ISET + 1 
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^ or>-> non o 


DO 40 J-1»NGRD 
00 20 KK«1» 100 
LINE (KK ) • 1H 
20 CONTINUE 

DO 30 KK«1,NGRD 
DO 30 LL "1>NGRD 
JND ■ NGRO* (KK - 1) ♦ IL 
KND ■ 10*(KK - 1) ♦ LL 
TF(KG( JND#J).NE.O) LINE (KND ) » 1HX 
30 CONTINUE 

IF(J.NE.l) WR IT E ( 6# 1 ) IZERO*LINE 
IF(J.EQ.l) WRITE ( 6# 2 ) IZEPO>LINE 
40 CONTINUE 

TSCAT IS THE TOTAL SCATTERING 
CROSS SECTION 

TSCAT ■ GLEG15(0.,PI,FUN) 

XNZA ■ 3.65E+20 

PTV - EXP(-TSCAT*XNZW*DIST) 

PTA - EXP(-TSCAT*XNZA*DIST) 

THE 170 LOOP IS OVER CONCENTRATIONS 

DO 170 I-1,NGRDP 

IF(I.NE.NGROP) TRANS ( I SET* I ) • 0. 

NHIT ■ 0 
N M I SS -0 
NSCAT • 0 
NNX » 0 
NNY - 0 
NN 7 ■ 0 
NNE ■ 0 

THE 150 LOOP IS OVER GAHHAS 

DO 150 K«1,NGAM 
THETAS ■ 0. 

PHIS - 0. 

THETAF • 0. 

PHI F • 0. 

FNFR • E7ER0 
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oooooooooo ooonnoon 


XI IS THE INITIAL LOCATION OF 
THE GAMMA ALONG THE X-AXIS* HERE 
ASSUMED TO BE UNIFORMLY DISTRIBUTED 
THE INITIAL Y-COORDINATE# YI* IS 
ASSUMED TO BE UNIFORMLY DISTRIBUTED 
ACROSS THF CONTAINER WALL 




XI ■ CDIAM+RANF (T ) 


THETAI IS THE INITIAL THETA- 
COMPONENT OF THE GAMMA VELOCITY* 
HERE ASSUMED TO BE UNIFORMLY 
DISTRIBUTED INSIDE THF 
COLLIMATION ANGLE. THE INITIAL 
PHI-COMPONENT IS ASSUMED TO BE 
UNIFORMLY DISTRIBUTED BETWEEN 
ZERO AND 2*PI 


C 

C 

c 

c 


c 

c 

c 

c 

c 

c 


THFTAI • -THETAL/2. ♦ 
PHI I • 2.*PI*RANF(T) 
THETA ■ THETAI 

PHI • PHI I 

ZI • 0. 


RANF (T )*THETAL 


YI ■ -CRAD + ?«*CRAD*RANF(T) 

XLEN ■ CDIAM/FLOATCNGRD) 

ZLEN - CLEN/ FLOAT (NGRD ) 

XI ■ XI + DIST*SIN(THETA)*CQS(PHI) 
YI » YI ♦ DlST*SIN{ THETA )*SIN (PHI) 
ZI ■ ZI + DIST*CQS (THETA ) 


IXI AND IZI ARE THE GRID INDICES 
OF THE CURRENT GAMMA LOCATION 

IXI • XI /XL FN ♦ 1, 

IZI - ZI/ZLEN + 1. 


THE FOLLOWING TESTS DETERMINE IF 
THE GAMMA HAS EXCEEDED THE X OR Y 
BOUNDARIES OF THE CONTAINER. 

FOR A CIRCULAR CONTAINER THEY 
COULD BE REPLACED WITH A TEST 
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c XI **2 + YI **2 • GT. CRAD**2 

C 

IF(IXI«6T.NGRD) GO TO 100 
TF(IXI.LT.l) GO TO 100 

IF{AB$(YI) .GT.CRAD) GO TO 110 
C 

C STOP IF ENERGY BELOW CUTOFF 

C 

IF(ENER.LE.EFINAL) GO TO 120 
C 

C STOP IF GAMMA OUTSIDE Z-BOUNDS 

C 

IF(IZI.LT.l) GO TO 130 
IF( IZI.GT.NGRD ) GO TO 140 
50 CONTINUE 
R ■ RANF(T) 

IXI - XI/XLEN + 1. 

IZI - ZI/ZLEN + 1. 

PT - PTA 

I F ( I • EQ • NGRDP ) GO TO 60 
JND » NGROMI-l) 

C 

C THE EFFECTS OF THE CONTAINER WALLS 

C COULD BE INTRODUCED HERE WITH THE 

C INTRODUCTION OF A »PT* FOR THE WALL 

C AS A FUNCTION OF XI, YI, AND ZI 

C 

IF(KG(JND+IXI»IZI).EQ.l) PT • PTW 
60 CONTINUE 
C 

C BRANCH IF NO SCATTER 

r 

IF(R.LT.PT) GO TO 90 
C 

C THE GAMMA SCATTERED 

r 

NSC AT • NSCAT + 1 
R - RANF(T) 

DO 70 KI » 1,179 

IF(YY{KI).LE.R.AND.YY(KI+1).GE.R) GO TO 80 
70 CONTINUF 
KI - 179 
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o n o r> o o o r> r> ooo 


c 

c 

c 


80 CONTINUE 


DETERMINE SCATTERING ANGLE 

THFTAS • XX(KI) + (XX(KH-l) - XX<KI>>*(R - YY (K I ) ) / 
1 ( YY ( K I + 1 ) - YY (KI ) ) 

PHIS - 2.*PI*RANF<T) 

THETA ■ THETA+180* /PI 
PHI • PHI+180./PI 
THETAS • THETAS+180./PI 
PHIS ■ PHI S+180 • /PI 

GET NEW DIRECTION 

CALL POST (THETA# PHI, THETAS# PHIS* THETAF,PHIF ) 

THETA • THETAF*P I / 180* 

PHI • PHIF*PI/180. 

THETAS • THETAS*PI/180. 

DECREMENT ENERGY 

ENER » ENER/C1. + ALPHA* ( 1 • - C0$( THETAS )> ) 

00 CONTINUE 


ADVANCF THE POSITION OF THE GAMMA 

XT ■ XI ♦ DIST*SIN(THETA)*COS(PHI> 

YT - YI + DIST*SIN(THETA)*SIN(PHI) 

ZI • ZI ♦ DIST*COSC THETA) 

IXI ■ XI/XLEN ♦ 1. 

IZI « ZI/ZLEN ♦ 1, 

REPEAT TESTS TO DROP THE 
GAMMA FROM FURTHER CONSIDERATION 

IF(IXI.GT.NGRD) GO TO 100 
IF(IXI.LT.l) GO TO 100 
IFtARS(YI).GT.CRAD) GO TO 110 
IF(ENFR.LE.EFINAL) GO TO 120 
IF(IZI.LT.l) GO TO 130 
I F ( I 71 • GT# NGRD ) GO TO 1A0 
GO TO 50 
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c 

c A MI SS> X OUTSIDE BOUNDS 

C 

100 NMISS ■ NMISS ♦ 1 
NNX • NNX + 1 

• GO TO 150 
C 

C A MISS* Y OUTSIDE BOUNDS 

* C 

110 NMISS - NMISS ♦ 1 
NNY ■ NNY + 1 
GO TO 150 
C 

C A MISS* FNERGY BELOW CUTOFF 

C 

120 NMISS - NMISS + 1 
NNE ■ NNE + 1 
GO TO 150 
C 

C A MISS* 1 OUTSIDE BOUNDS 

C 

130 Nmiss • NMISS + 1 
NNZ ■ NNZ + 1 
GO TO 150 
1 AO CONTINUE 
C 

C THE GAMMA REACHED THE FAR CONTAINER 

C WALL, HERE ASSUMED TO BE COUNTED 

C 

NHIT » NHIT + 1 
150 CONTINUE 
C 

C OUTPUT STATISTICS TO TAPE1 

C 

WRITE (1,3) I*NHIT*NMISS*NNX*NNY, NNZ*NNE 

IF(I.PO.NGRDP) GO TO 160 

TRANSf ISET,I) - FLOATt NHIT ) /FLOATtNGAM ) 

GO TO 170 

160 TTT(ISET) • F LOAT ( NHI T ) /FL OAT (NGAM ) 

170 CONTINUE 

WP I TF ( 6, A ) TTT( ISET), ( TR ANS ( IS ET# J ) , 1, NGRD ) 

1R0 RETURN 
END 
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SUBROUTINE DISP (TTT* TRANS* NS ET# NGRD ) 

OISP GRAPHICALLY DISPLAYS THE 
SUMMARY STATISTICS 

DIMENSION TTT(IOO) 

DIMENSION TRANS ( 100*10 ) 

DIMENSION X(15)*YU5) 

COMMON/RSLT/AVG(ll)*STDdl ) 

NPRDP - NGRD ♦ 1 
NGRDM ■ NGRD - 1 
X ( 1 > - 0. 

DO 10 I ■ 1* NGRD 

X < 1 + 1 ) - X C I > ♦ 100* /FLOAT (NGRD) 

10 CONTINUE 

DO 30 I » 2# NGRDP 
Y(I) • 0. 

DO 20 J-1,NSET 
Y { I ) - Y { I ) + TRANS ( J* I— 1 ) 

20 CONTINUE 

Y ( I ) - Y ( I ) /F LOAT (NSET) 

AVG(I) - Y(I) 

30 CONTTNUF 
Y ( 1 ) - 0. 

DO 40 J ■ 1* NSET 

Y ( 1 ) • Yd) ♦ TTT { J ) 

40 CONTINUE 

Y ( 1 ) • Y(1)/FL0AT(NSET) 

AVG < 1 ) ■ Y(l) 

X (NGRD+ 2 ) • 0. . 

X (NGRD+3) - 10. 

Y ( NGRD+2 ) » 0. 

Y (NGRD + 3) • 0.1 

CALL AXES(0.>0.,0.*10.*X(NGR0+2)»X( NGRD+3 ),1.,0.# 

1 7HT WATER*0.2*-7) 

CALL AXES(0..0.>90.>10.*Y(NGR0+2)> Y ( NGRD+3 )> l.*0.* 
1 12HTRANSMISSI0N>0.2#12) 

CALL LINPLT(X*Y>NGRDP*1#0>0#0*0) 

DO 50 J ■ l# NGRDM 
XP - 100+J /(10+NGRD) 

DO 50 I-1,NSET 
YP - TRANS(I» J)/0.1 
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CALL PNTPLT(XP#YP#21#1) 

50 CONTINUE 

DO 80 I-2,NGRDP 
STD ( I ) - 0. 

DO 60 J-1,NSET 

STD(I) ■ STD(I) + ( AV6 ( I ) - TRANS C J » 1-1 )) **2 
60 CONTINUE 

STD ( I ) - STD(I)/FLOATCNSET) 

STD ( I ) • SORT(STDCn) 

STD(l) - 0, 

DO 70 J ■ 1# NSET 

STD ( 1 ) • STDC1 ) + (AVG(l) - TTTCJ))**2 
70 CONTINUE 

STD(l) • STD(1)/FL0AT(NSET| 

STD ( 1 ) - SORT ( STD( 1) ) 

80 CONTINUE 
RETURN 
FND 
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SUBROUTINE POST <T1, PI# T2, P2# T3# P 3 ) 
DIMENSION A(4,B) 

CALL UROTOC(l,T2,A) 

CALL UAPPLY(0.*0.#1.#A,X # Y#Z) 

CALL UROTOC ( 3» P2> A ) 

CALL UAPPLY(X> Y,Z,A>X, Y, Z ) 

CALL UROTOC ( 1; Tl> A ) 

CALL UAPPLY(X,Y>Z,A,X,Y,Z) 

CALL UROTOC (3> Pl> A) 

CALL UAPPLY{X,Y,Z,A,X,Y,Z> 
P3»ATAN(X/Y)*1B0./3.14159265 
IF (X.LT.O..AND.Y.LT.O. ) P3-360.-P3 
IF ( X • L T • 0 • • AND • Y • GE .0 • ) P3»180»—P3 
IF (X.GE.O..AND.Y.GE.O.) P3«90.+P3 
IF (X.GE.0..AN0.Y.LT.0. ) P3—P3 
R»SQRT(X*X+Y*Y) 

IF (Y.LT.O.) R — R 
T3»ATAN(R/Z)*180./3.14159265 
IF (R.LT.O. .AND .Z.LT.O • ) T3-180.-T3 
IF (R.LT.O..AND.Z.GT.O.) T3«-T3 
IF (R.GE.O..AND.Z.GE.O.) T3-T3 
IF (R.GE.O. .AND. Z.LT.O. ) T3-180.+T3 
RETURN 
END 
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o o o r> 


SUBROUTINE UAPPLY (X#Y#Z# A# U#V#W) 

POST-MULTIPLIES A POINT VECTOR (X Y Z) BY 4*4 
TRANSFORMATION MATRIX A TO GIVE <U V W). 

RFAL A ( 4# 4 ) 

C H IS THE HOMOGENEOUS CO-ORDINATE# WHICH IS USED 
C TO EFFECT VARIOUS PROJECTIONS. IN THIS ROUTINE# 

C MATRIX A IS ASSUMED TO PRODUCE ONLY AFFINE 
C TRANSFORMATIONS# WHICH MEANS THE LAST COLUMN OF 
C MATRIX A IS (OOO 1)# SO H IS UNITY. BUT IT#S 
C HERE IN CASE ONE WISHES TO CHANGE IT... 

C (UNCOMMENT NEXT LINE AND DIVIDE X# Y# AND Z BY H.J 
C H « X*A( 1# 4) + Y*A ( 2# 4 ) + Z*A(3#4) «• A(4#4) 

U • ( X*A ( 1# 1 ) + Y*A (2,1) + Z*A ( 3# 1 ) + A(4,l)) 

V ■ ( X*A ( I# 2 ) + Y*A ( 2# 2 ) + Z*A(3#2) + A(4#2)) 

W « ( X*A { I# 3) «• Y*A( 2# 3 ) + Z*A(3#3) ♦ A(4#3)) 

RETURN 

FNP 



~>00000f>00 


SUBROUTINE UROTOC (IAXIS, ANGLE# A) 

CONSTRUCTS A 4*4 MATRIX THAT PERFORMS A ROTATION 
ABOUT X-, Y-, OR Z-AXIS OF ANY ANGLE. 

INPUTS — 

I AX IS 1-X 2«Y 3-Z 

ANGIE AMOUNT OF ROTATION, IN CURRENT UNITS. MEASURED 
POSITIVE BY RIGHT-HAND RULE. 

REAL A ( 4, 4 ) 

COMMON /TORAOS/ ANGFAC 
DATA ANGFAC / .01745329252/ 

IF ( I AX I S .LE. 0 .OR. I AXIS .GE. 4) GO TO 40 
CALL UCLR (A) 

A ( 4, 4 ) • 1. 

RAO » ANGFAC * ANGLE 
COSANG • COS (RAD) 

SINANG - SIN (RAD) 

GO TO (10, 20, 30), IAXIS 
10 A (1,1) . 1. 

A (2, 2 ) • COSANG 
A ( 3, 2 ) ■ - SINANG 
A ( 2, 3 ) » SINANG 
A (3, 3 ) • COSANG 
RETURN 

20 A ( 2, 2 ) • 1. 

A ( 1, 1 ) - COSANG 
A(3,l) • SINANG 
A(l,3) • - SINANG 
A ( 3, 3 ) - COSANG 
RETURN 

30 A ( 3, 3 ) - 1. 

A ( 1, 1 ) • COSANG 
A ( 2, 1 ) « - SINANG 
A (1, 2 ) - SINANG 
A ( 2, 2 ) ■ COSANG 
RETURN 

40 CALL UIDENT (A) 

PFTURN 

fnd 
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SUBROUTINE UCIR (A) 

C 

C CLEARS 4*4 MATRIX A TO ZEROS. 
C 

RFAL AI4,4) 

5 DO 10 I « 1, 4 

DO 10 J-1,4 
10 A(I,J) • 0. 

RFTURN 

FNO 


SUBROUTINE UIDENT (A) 

C 

C MAKES 4*4 IDENTITY MATRIX A. 

C 

PEAL A I 4 , 4 ) 

CALL UCLR (A) 

DO 10 1-1,4 
10 All, I) • 1. 

RETURN 

FND 


SUBROUTINE FUN I F,N ) 

DIMENSION F(l) 

COMMON P I, RO, AL PHA 
DP 10 I » 1, N 
THETA « F(I) 

ANS ■ 2.*PI*R0**2*SIN(THETA> 

TMP1 - 1. + ALP HA *11. - COSITHETA)) 
TMP2 « 1. + C0S(THETA)**2 
TMP3 - (1. - COSITHETA) )**2 
ANS • ANS*ll./TMPl)**2*ITMP2/2. ) 

1 *11. + ALPHA**2*TMP3/ITMP2*TMP1)) 
FID ■ ANS 
10 CONTINUE 
RFTURN 
FND 



o n r> o ,*> n n o o o o o o o n o o d n "l o o o r> 


FUNCTION GLEG15(A,B, FUNCT) 

FUNCTION GLFG15 EVALUATES INTEGRALS OF THE FORM 

B 

INTEGRAL CG(X)] OX 
A 

WHFRF A AND 8 ARE FINITE USING 15 POINT GUASS-LEGENOR E QUADRATURE. 
USE 

VALUF«GLEG15(A,B, FUNCT ) 

A THE LOVER LIMIT OF INTEGRATION. MUST BE FINITE. 

B THE UPPER LIMIT OF INTEGRATION, MUST BE FINITE. 

FUNCT. .THE NAME OF AN EXTERNAL SUBROUTINE TO CALCULATE THE 
VALUE OF G ( X ) AT SPECIFIED POINTS. 

SUBROUTINE FUNCT(F.N) 

F . .ON INPUT THE N VALUES OF X 
ON OUTPUT THE N VALUES G(X) 

THAT IS, FCn-GCFLIM I-1,...,N 
N..THE NUMBER OF FUNCTION EVALUATIONS. 
THE NAME GIVEN FUNCT MUST APPEAR IN AN EXTERNAL 
STATEMENT IN THE CALLING PROGRAM 
PIMFNSION W(8),R(8),F(15) 

DATA R/O. , .201194093997434, .394151347077563, 

* .570972172608538, .724417731360170, .848206583410427, 

* .937273392400705, .987992518020485/ 

DATA W/, 202578241925561, .198431485327111, .186161000015562, 

* . 166269 205 81699 3, .1395 70677926154, .10 7159220467171, 

* .703660474881081F-01, . 307532419961 172E-01/ 

C«( B-A)*.5 

D- (B + A)*.5 
F(1)-P 
DO 10 I-2,a 
S-P (I )*C 
F(i) -n+s 
10 F(i+ 7 )-n-s 

TALL FUNCT (F, 15) 


00000010 
00000020 
00000030 
00000040 
00000050 
00000060' 
00000070 
00000080. 
00000090 
00000100 
00000110 
00000120 
00000130 
00000140 
00000150 
0000C160 
00000170 
00000180 
00000190 
0 0000200 
00000210 
00000220 
00000230 
00000240 
00000250 
00000260 
00000270 
00000350 
00000360 
00000370 
00000360 
00000390 
00000400 
00000410 
00000420 
00000430 
00000440 
00000450 
00000460 
Q0G0047O 
00000480 
00000490 
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D"W( 1 )*F(1 ) 00000500 

00 20 I » 2> 8 00000510 

20 D"D+W(I)*(F(I )+F(I+7) ) 00000520 

GLE615-C+D 00000530 

PFTURN 00000540 

FND 00000550 


« 
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DISCUSSION 


As seen from Table III, the gamma ray transmissions for uniform planar source 
and the uniform linear source for various levels of fluid content in the container 
are essentially the same if the fluid distribution is assumed to be random. It 
might thus appear that a uniform linear source passing vertically through the 
center of the one side of the container would suffice for fluid content measure- 
ment. However, if the fluid is concentrated against one of the container walls or 
is confined to corner regions due to centrifugal forces, the linear source will 
not give a correct indication about the fluid content of the container. It is 
therefore concluded that a uniform planar Cs 137 source, coupled with an appropriate 
planar detector on the opposite face, would provide a more appropriate — though 
computationally less efficient basis for developing a universal nonintrusive monitor 
for liquid contents in sealed vessels. 

For the source/detector geometry illustrated in figure 1 and for a measurement 
accuracy of '< 10 percent, a planar Cs^- 37 gamma ray source of strength 10 pc will 
be adequate. Such an arrangement is expected to be particularly sensitive when 
the fluid content falls below 50 percent level. 


CONCLUDING REMARKS 

A nonintrusive nuclear monitor for liquid contents in sealed vessels has been 
described. Illustrative test studies indicate that the technique is capable of 
providing accurate (< +10%) information about the residual fluid content in sealed 
vessels. The monitor would be useable in all environments — ranging from 0-g of 
space to high-g of fast maneuvering aircraft. 
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Table I. Sample Formatted Display Output for Planar Source/Planar Detector Configuration 


Fluid 

content 


EZERO Ot 6620E+03 
EFINAL 0.2000E+03 
CLEN a 3048E+02 


CRAD 0. 1016E+02 
NSET 10 
NGRAD 10 


NGAM 400 

DIST a 1016E+01 

THETAL a 3000E+02 

XNZW a 3346E+24 


(%) 


Case 1 


Computed 

transmission-^ 


Case 2 


Computed 

transmission 


Case 3 


Computed 

transmission-* 


0 

10 

20 

30 

40 

50 

60 

70 

80 

90 

100 


X 

y yy 

y yy 

y yy 

y yy 

x yyy^ 

It x yyx 

x xx xxxxy 

X XXYXYXYX 

xxyyyyxxxx 



y 

y 

y y x 

yyy y x 

yyy x x y 

xxy y x y 

xxx xxx 

xxxxx yyy 

xxxxxxxxxx 


Y Y X 

y y xx 

y yyyy 

y y yyyy 

yyy yyxx 

xxx yyyy 

xxy yyyy 

yyy xxxx 

y y xxxxxy 

XXXXXXXXXX 


y 

y 

y y 

y yy 

y yyyy 

y yyyy 

x yyyx 

yxxyyxxxx 

XXYXXYXXXY 

xxxxxxxxxx 


y 

y 

y y 

yyy y 

yyy x 

y yxy y 

y xxx yyy 

xx xxy xxy 

xyxxxx xxx 

xxxxxxxxxx 


Y 

y 

y y 

y yy 

yyyy 

y y xyyy 

x yy yyyyy 

yyyy yxyxx 

xyxxxyyyyy 

xxxxxxxxxx 


y 

y yy 

y yy 

y yy 

y yyyyy 

y yyyxyy 

yyyyyxxyy 

XYYXyXXYY 

YyXYXXXXXX 

xxxxxxxxxx 


y 

Y 

y y 

y y 

x y 

yy yyy 

yyyy yyy 

YYXXXX X X 

xxxxxx xxx 

xxxxxxxxxx 



y 

y y y x y 

y vy yyy y 

y yyyyyy y 

y xyyyyyyy 

y yyyyyyxy 

y yyyyyyyy 

x xyyyyyyy 

XXXYXXXXXX 


y 

y y y y 

vyy x y 

yyyy y y 

xyxy y y 

yyyy x y 

yyyyy x v 

YXXXX X X 

XXXXX YXX 

YXYXXXXXXX 

.7*750" 

.607500 

• 5400CC 

.4*0000 

.407500 

♦350000 

.305000 

.225000 

#215000 

.170000 

.007500 


y 

y y y 

y y y y 

yy x y x 

yy x x x 

yx x yy x 

xxxx xx y 

yyxx xx x 

xxxxyxx x 

YXXXXXYXXX 


V 

y y 

yy yyy 

y vyv yyy 

y xyy yyy 

y xyy yyy 

y xxy xxx 

y xxyxxxy 

XYYXXYXXX 

YXXXXXYXXX 


y 

y 

y y 

yyy 

yy y y 

yyyyy yx 

xyyxyyyx 

y YXXXXYXy 

yyxxyyxxxx 

XXXXXXXXXX 

1 


y y 

yy y 

yy y x 

yy x x 

yy x y yy 

yy x yyyy 

yy x xxyxy 

xy y xxxxx 

xxxxxxxxxx 


y 

y 

y 

yyy 

yyy 

xxx x 

yyy xxx 

xxx xxx 

yyy yyyyx 

YXXYXXYXXX 


y y y 

yy y yy 

yy x xx 

yy x yy 

y yy y yyy 

X xx x xxx 

x yyyy yyy 

yyxyxx xxx 

xxxxxxxxxx 

xxxxxxxxxx 


yy 

yy v 

yy y 

y yy v 

y yy xx 

yx xx xx 

yy yyyyy 

yy xxyxxxy 

XX YYYXXXX 

xxxxxxxxxx 

; 

V 

y y 

X v 

y x 

xyy x 

yyx y 

Y xxx X 

x xxxx y 

XXXXYXX XX 

xxxxxxxxxx 



y 

y yy 

y y yy 

y y x y xx 

yx xxx x xx 

yyyyyyyyxy 

Yxxxxxvxxyi 

xxxyyyxYXx 

xxxxxxxxxx 




y x 

yyy 

yyy xy 

xxy yy 

xxx xx 

xxxxxyxx 

XYXXYXYXX 

xxxxxxxxxx 

.7f?*oo 1 

.*57*00 

.*17*00 

.417*500 

♦370000 

♦2O50CC 

.260000 

.167*100 

.100000 

.152500 

.110000 


y 

y y y 

y yy y 

yy xy yy 

yy xx xx 

yx x xx xy 

yy y xy xy 

yy x xxxxx 

yy yyyyyxy 

xxxyxxxxxx 


y 

yy 

y yy 

y y yy 

yxy x xyy 

yyy yyxyy 

yyy yyyyyy 

yyy xxxyyy 

XXX yXYYXY 

xxxxxxxxxx 



y 

yy y 

yyyy y 

yxyyx y 

yyyyy yy 

yyyyy xxx 

XXXXX YYXXi 

XXXXXXXXXX 

xxxxxxxxxx 


y 

y 

y y 

x y 

y y 

y y x 

xxy xx 

yyy xx xx 

yxxxyy xx 

xxxxxxxxxx 


y 

y y 

y y 

yy y 

yy x yy 

yyyyy xx 

XXXXYYXX 

y yxyxxxxx 

X XXXYYXXX 

xxxxxxxxxx 


Y 

y 

y 

y yy 

y yyy 

yyy yyxx 

xxxx yxxx 

YXXXXYXYXX 

yyyxxx xxxx 

xxxxxxxxxx 




yy y y 

yy > y 

yy yyx 

xx xxx 

xx y x y x 

yy x x y x 

XYXXXYXX X 

xxxxxxxxxx 


y 

y y y 

y Y YY 

y y xx 

x yy y xx 

v yy x xx 

x xxyyy xy 

yxxxxxx xx 

XYXXXYXXXY 

xxxxxxxxxx 


X 

y x 

y x 

yyy 

yyy 

x . x x 

y x xy 

xx xx xx 

yyy xy xxx 

xxxxxxxxxx 


y y y 

yy y yy 

yy y yy 

yy y yyy 

yy y xvy 

yy y yyxx 

xx y xxxy 

XX X YXYX 

yy xyyxyx 

xxxxxxxxxx 

* 75? 500 

#66?600 

.560000 

.475000 

#3*0000 

.317500 

• 305 OCC 

.237500 

.185000 

.157500 

.130000 




Table II. Sample Formatted Summary Output for Planar Source/Detector 
Configuration. 

EZERO 0.6620E+03 
EFINAL 0. 2000E+03 
CLEN 0.3048E+02 

CRAD 0.1016E+02 

NSET 10 

Gamma Ray 

Pe rcent Water Transmission Std. Dev. 


NGRAD 

10 

NGAM 

400 

DIST 

0. 1016E+01 

THETAL 

0. 3000E+02 

XNZW 

0. 3346E+24 


0.0000 

0.7635 

0.0221 

10.0000 

0.6420 

0.0324 

20.0000 

0.5350 

0.0317 

30.0000 

0.4538 

0.0264 

40.0000 

0.3795 

0.0220 

50.0000 

0.3090 

0.0239 

60.0000 

0.2705 

0.0302 

70.0000 

0.2263 

0.0228 

80.0000 

0.1908 

0.0177 

90.0000 

0.1583 

0.0107 

100.0000 

0.1225 

0.0209 



Table III. Summary of Transmitted Gamma Ray Intensities 
for Various Fluid Contents. 

(a) Uniform Planar Source 


Percent Water 

! 

; Transmission 

(*) 
R V 1 

(%) 

j Random 

; Fluid Distribution 

] . - - 

Non-Random 
Fluid Distribution 

Random 

Fluid Distribution 

Non-Random 
Fluid Distribution 

0.0000 

• 0.7635+0.0221 

0.7560+0.0198 

6.233+0.200 

5.458+0.133 

10.0000 

0.6420+0.0324 

0.6002+0.0175 

5.241+0.221 

4.334+0.136 

20.0000 

• 0.5350+0.0317 

1 

0.4777+0.0193 

4.367+0.230 

3.449+0.147 

30.0000 

; 0.4538+0.0264 

0.3890+0.0271 

3.704+0.229 

2.809+0.177 

40.0000 

0.3795+0.0220 

0.3103+0.0140 

3.098+0.229 

2.240+0.152 

50.0000 

! 0.3090+0.0239 

0.2670+0.0131 

2.522+0.248 

1.928+0.156 

60.0000 

i 

; 0.2705+0.0302 

0.2343+0.0205 

2.208+0.282 

1.692+0.194 

70.0000 

: 0.2263+0.0228 

0.2025+0.0165 

1.347+0.271 

1.462+0.183 

80.0000 

; 0.1908+0.0177 

0.1755+0.0149 

1.558+0.263 

1.267+0.192 

90.0000 

| 0.1583+0.0107 

0.1465+0.0159 

1.292+0.238 

1.058+0.215 

100.0000 

0.1228+0.0209 

0.1385+0.0148 

1 . 000+0 . 341 

1.000+0.214 


1 

i 

1 

i — 

i 

1 


(*) _ Transmission for x % Water 

Transmission for 100% Water 



Table III. Summary of Transmitted Gamma Ray Intensities 
for Various Fluid Contents. 

(b) Uniform Linear Source (*) 


Percent Water 

Transmission 


R^ 

*) 

(%) 

Random 

Fluid Distribution 

Non-Random 
Fluid Distribution 

Random 

Fluid Distribution 

Non-Random 
Fluid Distribution 

0.0000 

0.8732+0.0201 

0.8682+0.0134 

5.940+0.158 


5.788+0.115 

10.0000 

0.7407+0.0289 

0.7080+0.0184 

5.039+0.174 


4.720+0.126 

20.0000 

0.6327+0.0281 

0.5727+0.0153 

4.304+0.179 


3.818+0.127 

30.0000 

0.5277+0.0303 

0.4560+0.0243 

3.590+0.192 


3.040+0.153 

40.0000 

0.4455+0.0246 

0.3653+0.0203 

3.031+0.190 


2.435+0.156 

50.0000 

0.3770+0.0320 

— — — — 

0.3070+0.0177 

2.565+0.220 


2.047+0.158 

60.0000 

0.3295+0.0182 

0.2660+0.0174 

2.241+0.190 


1.773+0.165 

70.0000 

0.2613+0.0196 

0.2293+0.0227 

1.778+0.210 


1.529+0.159 

80.0000 

0.2088+0.0280 

0.1955+0.0152 

1.420+0.269 


1.303+0.178 

90.0000 

0.1785+0.0177 

0.1640+0.0117 

1.214+0.234 


1.093+0.171 

100.0000 

0.1470+0.0198 

— 

0.1500+0.0150 

1.000+0.269 


1.000+0.200 


(*) R Transmission for x % Water 
Transmission for 100% Water 



FIGURE L - SCHEMATIC DIAGRAM FOR THE UNIFORM 

PLANAR SOURCE/PLANAR DETECTOR CONFIGURATION. 
THE CONTAINER IS ASSUMED TO BE A lx 2d x 2d 
RECTANGULAR BOX FOR COMPUTATIONAL PURPOSES. 
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